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Metastability of a particle trapped in a well with a time-periodically oscillating barrier is studied 
in the Floquet formalism. It is shown that the oscillating barrier causes the system to decay faster 
in general. However, avoided crossings of metastable states can occur with the less stable states 
crossing over to the more stable ones. If in the static well there exists a bound state, then it is 
possible to stabilize a metastable state by adiabatically increasing the oscillating frequency of the 
. barrier so that the unstable state will eventually cross over to the stable bound state. It is also 

found that increasing the amplitude of the oscillating field may change a direct crossing of states 
f*"^ ' into an avoided one. 



PACS numbers: 03.65.Xp, 33.80.Be, 74.50.+r 



INTRODUCTION 



Ever since the advent of quantum mechanics, quantum tunneling has been an important and fascinating subject. 
This phenomenon arises frequently in physics. In fact, one of the first successful applications of quantum mechanics has 
been the explanation of the a decay of atoms as a quantum tunneling process pj . Recent examples include tunneling 
phenomena in semiconductors and superconductors Q, in Josephson junction systems 0, resonant tunneling in 
heterojunction nanostructures Q, tunneling ionization of atoms p|, photon-assisted tunneling in superconducting 
•^-j- ' junctions and semiconductor superlattices [6(, etc. 

In cosmology, quantum metastable states play an essential role in some versions of the inflationary models of the 
■ early universe y\ . In these models inflation of the early universe is governed by a Higgs field trapped in a metastable 
state. Inflation ends when the metastable state decays to the true ground state of the universe. During inflation 
r-j 1 the universe expands exponentially. It is thus obvious that the metastable state of the Higgs field is trapped in a 
O ,' rapidly varying potential. The problem is therefore a truly time-dependent one. Unfortunately, owing to the inherent 
difficulties of the problem, more often than not one has to consider the decay of the Higgs field in a quasistationary 
approximation, in which the decay is studied assuming a static potential |8|]. Surely this approximation is hard to 
justify, but for the present one has to be content with it. Ultimately one hopes to be able to tackle the nonstationary 
^ ■ case. To this end, it is desirable to gain some insights first by studying metastability in time-dependent potential in 
| simple quantum-mechanical models. 

An early attempt at studying the effects of time- varying forces on quantum metastability appears in Fisher's work 
k>( \ 13 j which was motivated by an experiment on quantum tunneling of the phase in a current-biased Josephson junction 
with a weak microwave perturbation |Toj . In this work Fisher considered the general problem of quantum tunneling 
in a metastable well with a weak oscillatory force. There he reformulated the standard WKBJ approach to quantum 
decay in order to include a weak time-dependent perturbation. For a class of metastable potentials which interpolates 
between the cubic potential and a truncated harmonic-oscillator potential, he showed that the decay rate is generally 
enhanced by the weak oscillatory force. The potential considered by Fisher has a number of oscillator-like levels near 
its minimum. The opposite situation where only two levels are present was considered by Sokolovski [Tl| , who studied 
the effect of a small ac field mixing two levels in the well on the tunneling rate in a scmiclassical framework. 

The results in 0, are quite general for a class of weak oscillatory forces. However, it is desirable to consider 
other possibilities, e.g., exact solutions and/or nonperturbative results. In this respect, we had considered previously 
M an exactly solvable quantum metastable system with a moving potential which has height and width scaled in a 
specific way introduced by Berry and Klein [jjj . We found that in this model a small but finite nondecay probability 
could persist at large time limit for an expanding potential. 

In this paper we consider another simple driven quantum metastable model in which a particle is trapped in a well 
with a periodically driven rectangular barrier. In order to do away with any restriction of amplitude or frequency of 
the driving force, and of the number of states in the potential, we treat the problem in the framework of the Floquet 
formalism [lj, fl5L \WL Il7| . An exact expression determining the Floquet quasienergies of stable or metastable states 
in the well is derived. From the solution of this equation we find that while the oscillating barrier makes the system 
decay faster in general, there is the possibility that avoided crossings of metastable states can occur with the less 
stable states crossing over to the more stable ones. 
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That an oscillating potential can affect the tunneling property of a system has also been noticed before, e.g., in 
quantum transport process 0, 0, 0, . In 0, 0] it was found that a particle can be localized in one side of a 
time-dependent double well if the amplitude and the frequency of the oscillating field were chosen properly. In it 
was demonstrated that a propagating particle at appropriate incident energy can be trapped into a bound state by 
an oscillating square well. Our example shows how a time-periodic field can modify the metastability of a decaying 
state. 



II. THE MODEL 




The model we consider consists of a particle of mass m trapped in a square well with a harmonically oscillating 
barrier, 

x < , 

1 ! ' ' ' 1 = ^ vl j. Vi cos(wi) , a < x < b ,' ^ 

x > b . 

Here Vo, V\, Vq and uj are positive parameters, with V ' < Vq and V\ < Vo — Vq. According to the Floquet theorem, the 
wave function of a time-periodic system has the form ^ e (x, t) — e~' let / h Q E (x, t), where & e {x, t) is a periodic function 
with the period T = 2tc/lj, i.e., & € (x, t + T) — $ s (x, t), and e is the Floquet quasienergy, which we will call Floquet 
energy for brevity. It should be noted that the Floquet energy is determined only modulo %uj. For if {e, $ E } is a 
solution of the Schrodinger equation, then {e' = e + nhuj , $ e < = $ £ exp(zno;t )| i s also a solution for any integer n. 
But they are physically equivalent as the total wave function \l/ e is the same [lq . All physically inequivalent states 
can be characterized by their reduced Floquet energies in a zone with a width Tiuj. We therefore consider solutions of 
e only in the first Floquet zone, i.e., e € [0, h u>). 

Following the procedures described in [t(| (see also 0|), we get the wave function as follows: 

*(a?,t) = e- iet / h $ £ {x,t) (2) 

E" = -oo A n sm(k n x)e- m ^ , < x < a , 

EZ-oo EZ-oo (W q,x + bie- qix ) Jn-l (Vx/hco) e~ m ^ , a < x < b , 

EZ-oo t n e< x e- inut , x > b 



■iet/h 



where 



k n = \/2m(e + nhuj)/Ti , 

qi = ^2m(V -e- lhw)/h , (3) 
k' n = yj 2m(e + nfku - Vq)/U , 

and J n is the Bessel function. In the region x > b, we have adopted Gamow's outgoing boundary condition, namely, 
there is no particle approaching the barrier from the right 0]. Matching the wave function and its first derivative at 
the boundaries x — a and x — b, we obtain the relations among the coefficients A n , a n , b n and t n : 

A n sm(k n a) = Y,( a i eq ' a + b i e ~ q ' a ) Jn-i{a) , 
i 

k n A n cos(k n a) = {aie qia - M^ ) J n -i{a) , (4) 

i 

t n e lk '" b - Y,( aie9 ' b + bie ~ q ' b ) J ^ a ) ' 
i 

where a = Vi/Tuo. The Floquet energy is determined from these relations by demanding nontrivial solutions of the 
coefficients. In practice, however, we must truncate the above equations to a finite number of terms, or sidebands as 
they are usually called in the literature, eg., n = 0, ±1, . . . , ±JV . The number is determined by the frequency and 
the strength of the oscillation as N > V\/hu fig 1 . 
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We proceed to determine the Floquet energy as follows. We first separate the boundary conditions for the central 
band (n = 0) from those for the subbands (n ^ 0) in Eq.QJ. From the boundary conditions for the subbands (n ^ 0), 
one can relate the coefficients ai and bi (I / 0) with the coefficients do and bo through the following two relations: 

ai = fia(ko, k' , lu, V\)a a + fw(ko, k' , to, Vi)b , (5) 
bi = gi a (k ,k' ,u;,Vi)ao + gib(k ,k' ,u;,Vi)bo , (6) 

where /'s and g's are functions determined as follows. Eliminating the A n 's and t n 's in Eq.(0}, we can obtain 
A- n e^ a J Q a n + Al n e-^ a J b n + £ A~^ a J n „ m + £ A+^J^k 

= -A- Q ei° a J n a - A+ e~ q ° a J n b , (7) 

and 

B+ n e^ b J a n + B- n e-^ b J b n + £ B+ l e^ b J n . l ai + ^ B^e^J^k 
= -B+ e q " b J n a - B- e- q ° b J n b , (8) 

with 

A^ t = cos k n a ± j- sin k n a , and B^ t = 1 ± z-^- . (9) 

Equations J7J) and (JSJ) allow us to solve for a; and bi in terms of do and bo in the forms of Eqs.© and © by means 
of the Cramer's rule in matrix algebra. As mentioned before, in practice a truncated version of Eqs. (JTJ and (|SJ has 
to be used. 

Using Eqs.JHJ an( l © we can rewrite the boundary conditions for the central band n = as 

A sm(k a) = F 1 (k ,k' ;u;,V 1 )e^ a a + F 2 (ko,k' ;iu,V 1 )e-^ a bo , 
koA cos(k a) = F^ko^^u^V^qoe^ao-Fiiko^^Lu^qoe-^bo , (fO) 
t e lk '° b = F 5 (koX;u,V 1 )e« ob a + F 6 (koX;u,V 1 )e-' lob b , 
ik' t e< b = F 7 (k ,k , ;u;,V 1 )qoe llab a ^F a (k ,k' ;iu,V 1 )qoe' qob bo , 
where the coefficients Fi(k ; u 7 V\) (i = 1, . . . , 8) are 

F 1 (k X;«>,Vi) = J (a) + e-v° a J2(fiae qia +g la e-« a )j^(a), 

F 2 (k ,k ;u,V 1 ) = Ma) + e q ° a Y,{fibe qia +9ibe- qia )J-i{a) , 

F 3 (fcoX;w,^i) - Jo(a) + e^ oa V-(/ ; ae 9iQ -ff ;a e- 9ia ) J_,(a) , 

F 4 (fc ,/^,lq) = J (a)-e«°y)^.(/,6e* -<R 6 e-* ) J_j(a) , 

Fsfe^iw,^) = Jo(a)+e-* 6 5^(/ Io c« 6 + fl Ia e-* 6 ) J_,(a) , 

f«(fcoX;w,^i) = Jo(a)+e 9ob ^(/ ib e* b +. 9ib e-* b ) J_,(a) , 

F 7 (fco,^;w,^) = Jo(a)+ e -* b V * (/ iQ e* b - 5iae -«' b ) J_ ; («) , 

FaikoXwVi) = J Q (a) - e qob T^ (f lb e qib - g lb e- q ' b ) J_ t (a) . 

By demanding nontrivial solutions of the coefficients ao, bo, Aq, and to in Ea. (|f Ojl . we obtain an equation which 
determines the Floquet energy e: 



F4—t&nk a + F 2 = — — — [ F 3 — tanfc a-Fi e 9lH J 



(11) 
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We recall here that ko,qo,k' are functions of the Floquet energy e [c.f. Eq. J5J]. If the solutions e of Ea jllj l are 
complex (real) numbers, the corresponding Floquet states are metastable (stable) states. The nondecay probability 
P(t), which is the probability of the particle still being trapped by the potential barrier at time t > 0, is given by 



P(t) 



tf\*(x,0)\ 2 dx 

= c 2Im(e)t/h fo \* S (x,t)fdx 

/ b |$ £ (x,0)| 2 dx 

_ e 2Irn(e)t/h h ^ ^ 

with P(0) = 1. The imaginary part of the Floquet energy, which enters P(t) via the factor exp[2Im(e)t/Ti], gives a 
measure of the stability of the system. Unlike the static case, however, here P(t) is not a monotonic function of time, 
owing to the time-dependent function h(t) after the exponential factor in Ea. H12|) . But since h(t) is only a periodic 
function oscillating between two values which are of order one, the essential behavior of P(t) at large times is still 
mainly governed by the exponential factor. Hence, as a useful measure of the nondecay rate of the particle in the 
well, we propose a coarse-grained nondecay probability P(t) defined as 

P(t) = e 2Im{£)t/h (h(t)) , (13) 

where (h(t)) is the time average of h(t) over one period of oscillation. The graphs of P(t) and P(t) will be given in 
the next section. 

It is easily seen that the coefficients Fi(ko;cu, V\) all approach one in the limit a — Vi/Tico — > 0, 

lim Fiiko-^^Vx) — > 1 , i= 1,...,8 . (14) 

Vx/u>— >0 

Hence in the limit V\ — > or uj — > oo, Eq.JTTJ reduces to the corresponding equation for the case of static potential with 
potential Vb in the region a < x < b, and the Floquet energy in this limit is just the (real or complex) eigenenergy of 
the static case. This is understandable, since in the limit V\ — > the potential becomes static, and at high frequencies 
the particle in the well will only see a time-averaged barrier of effective height Vb 0] . 



III. NUMERICAL RESULTS 



We now study numerical solutions of Ea. (|llf) with a specific potential. We take a = 1, 6 = 2, Vb = 15, and 
Vq = Vb/2 in the atomic units (a.u.) (e = m e = h = 1). In the static case this potential supports one bound state, 
with energy E /Vq = 0.232123, and one metastable state, with complex energy Ei/Vo = 0.864945 - 0.0025526U. For 
the oscillating potential, we solve Ea. in 2-sideband approximation, i.e., we take N = 2. This is accurate enough 
for oscillating frequency u) > Vi/2. 

Figures 1 and 2 present the graphs of the real and imaginary parts of the Floquet energy (e/Vb) as a function of 
oj/Vq > 0.2 with V\ — O.lVb and 0.2Vb, respectively. We find that the solutions of Ea. i|ll|) have the form e — Eq + nu 
(n = 0, ±1, ±2, . . .), with Re(eo) (the horizontal branch) laying close to the energies Eq and Re{E\) in the static 
potential. That is, these branches of Re(e) emanate from either Eq or Re(Ei) at u> = 0. Branches emerging from 
the same point have the same imaginary part. Numerical results show that, with the barrier oscillating, the stable 
state (Eq) in the static case becomes unstable, and the unstable state (Ei) will decay even faster. For simplicity, in 
Figs, (la) and (2a) we show only six branches (n = 0, ±1, ±2 and —3) emerging from Re(Ei), and only the central 
branch (n — 0) and a subband (n = —1) from Eq. As mentioned before, we only take solutions in the first Floquet 
zone, Re(e) (modulo u), which are points under the line Re(e) = u>. 

In Fig. 3 we give the graphs of the probability density | \3/ 1 2 in the well and the barrier with the same parameters 
as in Fig. 2 for the two metastable states at frequency lo/Vq = 0.62. The Floquet energies of the less stable and the 
more stable state are e/Vq = 0.251714 — 0.004995i and 0.227343 — 0.001456i, respectively. Four time frames, namely, 
t = t X Vq = 0, 100, 200 and 300 a.u., are shown, with the probability density normalized to unity within < x < 2 
at r = 0. One can clearly see that the less stable state (dashed curve) decays much faster than the more stable state 
(solid curve). The nondecay probability (oscillatory curve) P(t), Ea. l|12|l . and the coarse-grained nondecay probability 
(monotonic curve) P(t), Eq. (|13fl . of these two states are shown in Fig. 4. It is clear that the coarse-grained function 
P(t) is a monotonic function, and does give a smooth measure of the stability of the system. 

From Figs. 1 and 2 we also see that a direct crossing occurs at frequency lo k, Re(E\ — Eq)/2 (point c). However, as 
u> approaches the frequency u> w Re(E\ — Eq) = 0.632822Vb, an avoided crossing (e, e') between the real parts of the 
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Floquet energies occurs. Figure 2 indicates that larger values of V± only enhance the instability of the system and the 
repulsion between the two levels at avoided crossing. Thus as the frequency u> is increased, the state emanating from 
Ei has Floquet energy with real part given by values along the path abb'cdd'e'f (the dark dotted curve), while the 
real part of Floquet energy of the state emerging from Eq lies along the path cegg'h (the solid curve). The imaginary 
parts of these two paths are depicted in Figs, (lb) and (2b). One sees that an exchange of the imaginary parts takes 
place at the avoided crossing ee'. In Figs. 5 and 6, we show the probability density of the two states within the 
potential barrier just before and just after the avoided crossing. Together with Fig. 3, these plots demonstrate clearly 
the switching of the states. Beyond the avoided crossing, the upper state becomes more stable than the lower state. 
This gives the possibility of stabilizing an unstable state by an oscillating field. We recall that as u> — > oo, Eg. (|11|> 
reduces to the one for the static potential. In the example considered here, the lower state supported by the static 
well is a stable bound state, and hence the unstable upper Floquet state can be made stable in the high frequency 
limit. Even more simply, the same aim can be achieved by adiabatically tuning down the amplitude Vi just after the 
avoided crossing, as in this limit the potential becomes the static one. 

Finally, it is interesting to note here the role of amplitude Vi of the oscillating barrier in the model. As we have 
seen, the presence of a nonvanishing Vi always makes the system less stable. However, if Vi is reduced, an avoided 
crossing may turn into a direct one. In the present case, the avoided crossing ee' changes into a direct crossing for 
Vi/Vq < 0.03. Conversely, increasing Vi could change a direct crossing into an avoided crossing. At an avoided 
crossing, the imaginary parts of the Floquet energies cross, while the real parts do not. At a direct crossing, it is the 
real parts, not the imaginary parts, that cross. But the more (less) stable state has the tendency to become less (more) 
stable. This is evident from the Floquet energy at the direct crossing point c in Figs. 1 and 2. These observations 
are consistent with the semiclassical results obtained in |l l| by perturbative methods. Hence, by a combination of 
adiabatic changes of the frequency and the amplitude of the oscillating barrier, one can manipulate the stability of 
different states in a quantum potential: tune up Vi until a direct crossing becomes an avoided one, increase u> so that 
the avoided crossing is passed, then reduce Vi to make the potential static. In the process, two states in the well are 
interchanged. 

IV. SUMMARY 

To summarize, our results show that an oscillating potential barrier generally makes a metastable system decay 
faster. However, the existence of avoided crossings of metastable states can switch a less stable state to a more stable 
one. If in the static well there exists a bound state, then it is possible to stabilize a metastable state by adiabatically 
changing the oscillating frequency and amplitude of the barrier so that the unstable state will eventually cross over 
to the stable bound state. Thus a time-dependent potential can be used to control the stability of a particle trapped 
in a well. 

Finally, we would like to comment on the differences between the stabilization of the decaying state discussed 
in this work and the interesting phenomenon of the suppression of ionization of atom (also called stabilization of 
atom) in superintense, high frequency laser pulses [20|. In this later phenomenon, it was found that while an atom 
is generally ionized by absorption of photons, the ionization rate of the atom can be dramatically suppressed as the 
intensity (amplitude) of the laser field exceeds a certain threshold. Thus stabilization of the atom is attained by 
increasing the intensity. On the contrary, stabilization of the decaying state described here is achieved by increasing 
the frequency of the oscillating barrier to the threshold at which an avoided crossing in the Floquet energies takes 
place, regardless of the amplitude of the field. Also, in our case stabilization is against quantum tunneling through 
a potential barrier, while in the case of atomic stabilization it is against ionization by photon absorption. Another 
difference of the two phenomena is that an atom in the ground state is stable when not being irradiated by the laser 
field, but the decaying state to be stabilized in our system is already unstable even in the absence of the oscillating 
field owing to quantum tunneling effect. Furthermore, suppression of ionization can be studied using the methods of 
classical nonlinear dynamics and chaos 21], but tunneling through a barrier considered here is a characteristic and 
fundamental quantum phenomenon. 
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Figure Captions 

Figure 1. The Floquet energies (e/Vo) of the two metastable states versus the barrier oscillating frequency (uj/Vq) 
for Vo = lba.u.,V^ = Vq/2 and V\/Vq = 0.1 in the atomic units (a.u.) (e = m e = U = 1) . In (a) the 
real parts of the Floquet energies are shown in the first Floquet zone under the line Re(s) = uj (the straight 
line). The light dotted lines show how the different branches emanate from the two states in the static case 
[with Eq/Vq = 0.232123 and Re(E 1 )/V = 0.864945]. In (b) the corresponding imaginary parts of the Floquet 
energies of the two states are plotted. The dotted curve corresponds to the state with real parts given along the 
path abb' cddl e' / , and the solid curve corresponds to the state with real parts given along cegg'h. 

Figure 2. Same plot as Fig. 1 for V = lba.u., Vq = V&/2, and V 1 /V a = 0.2. 

Figure 3. Probability density l^j 2 in the well and the barrier with the same parameters as in Fig. 2 for the two 
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metastable states at frequency lu/Vq = 0.62. Probability density normalized to unity within < x < 2 at 
t = Wq = a.u.. The more (less) stable state is indicated by solid (dashed) curve. 

Figure 4. Nondecay probability (oscillatory curve) P(t) and the coarse-grained nondecay probability (monotonic 
curve) P(t) as a function of time for the more stable (a) and the less stable (b) state in Fig. 3. 

Figure 5. Same plot as Fig. 3 at frequency lo/Vq — 0.63, just before the avoided crossing. The lifetimes of these two 
metastable states are comparable. 

Figure 6. Same plot as Fig. 3 at frequency oj/Vo = 0.64, just after the avoided crossing. The dashed (solid) curve 
represents the originally less (more) stable state, which now becomes more (less) stable. 
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Figure 5 
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